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We report detailed measurements of the pump-current dependency of the 
Oh 



self-pulsating frequency of semiconductor CD lasers. A distinct kink in this 



dependence is found and explained using rate-equation model. The kink 
denotes a transition between a region where the self-pulsations are weakly 
sustained relaxation oscillations and a region where Q-switching takes place. 
Simulations show that spontaneous emission noise plays a crucial role for the 
cross-over. 



I. INTRODUCTION 

Self-pulsating semiconductor lasers (SPSL's) are of great interest owing to their potential 
application in telecommunication systems as well as in optical data storage applications. In 
particular, in the latter case they are realized as so-called narrow-stripe geometry CD lasers 
where the self-pulsation is achieved via saturable absorption in the transverse dimension 
limiting the active region. A profound knowledge and understanding of their operation 
dynamics is therefore desired. 

SPSL's have been studied since the first diode lasers became available in the late 1960s 
|J. These first semiconductor lasers, although designed to operate in continuous wave 
(CW) mode, showed self-induced pulsations of the light intensity due to a combination of 
two reasons: (i) the laser resonance is internally excited through the nonlinear interaction 
of various longitudinal laser modes, thus causing mode beating at very high frequency; 
(ii) defects in the active material act as saturable absorbing areas, thus causing absorptive 
Q-switching processes. 

In the case of self-pulsations caused by saturable-absorbing effects, the self-pulsation 
frequency (SPF) dependence on the pump current was investigated in 0. In later works the 
self pulsations were attributed to undamped relaxation oscillations (RO) [@,f|. The precise 
values of the ROF, as calculated from a small-signal analysis, and the actual SPF, highly 
nonlinear, are however different, the SPF being always smaller than the ROF ||. 

Saturable absorption effects, causing self-pulsations in stripe-geometry lasers have been 
investigated since the early 1980s ||-f§. Saturable absorption is also responsible for self- 
pulsations in double-section laser diodes ||. A similar mechanism of dispersive Q-switching 
has been invoked to describe self-pulsations in multisection Distributed Feedback Lasers 



In this paper we study both experimentally and theoretically the dependence of the self- 
pulsation frequency (SPF) of narrow-stripe geometry self-pulsating semiconductor lasers, 
also known as CD-lasers, on the bias pump current. In these lasers, self-pulsation is induced 



via saturable absorption in the transverse dimension of the active region. The rate-equation 
model of Ref. H has been proven to be quite successful in describing the mechanism of 
self-pulsation and has already been used with success in analyzing such lasers subject to 
weak optical feedback ||. There it was found that, with and without feedback, there are 
two distinct regions in the SPF vs. pump-current curve, one where spontaneous emission 
dominates the laser dynamics between pulses and one where spontaneous emission always 
plays a minor role. 

In section II we present detailed measurements of the SPF vs pump-current curve. This 
curve confirms most of the findings of ||, and also shows a distinct cross-over point dis- 
tinguishing between linear and square-root-like behavior. In section III we confront the 
experimental results with a theoretical model, inspired by Ref. . Its results agree qualita- 
tively well with the experimental results, showing a distinct cross-over region. The location 
of the cross-over region is shown to be determined by the spontaneous emission rate. In 
Section IV we discuss the relationship between the SPF and the ROF using a small signal 
analysis. We discuss the various bifurcations that are predicted by our model, and compare 
it with the model of Ref. |J . 

II. EXPERIMENT 

We use a SHARP CD semiconductor laser diode, model LT022MD. The laser emits a 
continuous train of regular pulses with a frequency that depends on the bias pump current. 
A bulk layer of AlGaAs constitutes the active layer of this Fabry-Perot cavity that emits 
at 800-nm wavelength. The gain section is defined by the p-electrical contact and has the 
following approximate dimensions: 250-/xm long, 2-/im wide, and 0.2-//m thick. A very 
narrow contact of ~ 2/xm allows for current injection. Since the region capable of stimulated 
emission extends to both sides beyond the narrow stripe of the current contact, the wings 
of the optical field distribution will interact with these unpumped, and therefore absorbing, 
regions. In fact, these regions are saturably absorbing; when the optical intensity in the 



wings of the mode is large enough, the electron-hole pair population in the unpumped region 
reaches transparency, thus allowing a "self-Q-switched" pulse. There is no sharp boundary 
between the pumped and unpumped regions, making carrier diffusion an important effect. 
Indeed, in the model of Ref. || carrier diffusion between the pumped and unpumped regions 
is crucial for the appearance of self-pulsation. 

The experimental set-up is illustrated in Fig. 1. The laser is temperature controlled by 
a Peltier cooler at 20 C. It is DC biased by a low noise current supply. Laser emission is 
collected by an anti-reflection coated 0.65 N.A. laser diode lens. The resulting parallel beam 
is passed through a 30-dB isolator to avoid spurious effects caused by optical feedback and 
is launched into a 60-GHz photodiode (New Focus Model 1006). The converted electrical 
signal is observed with a 22-GHz bandwidth spectrum analyzer (HP 8563A). The typical RF 
spectrum of the SPSL is characterized by a main peak at the SPF, followed by overtones. The 
uncertainty on the self-pulsation frequency measurement is mainly due to the measurement 
of bias current, which has an error of less than O.f mA. The resolution of the spectrum 
analyzer is 100 kHz and the video filter is 30 kHz. 

For values near the threshold current, the low power emission makes it difficult to observe 
the signal. The value of the spectral density of the self pulsations is very close to the noise 
level and also the width of the feature in the power spectrum is wider than at higher currents. 
To overcome this problem, a small current modulation is applied to the device for injection 



currents below 47 mA |[13| , [[LJ]. Its power is kept sufficiently low so that it does not affect 
the oscillation behaviour of the laser and does not induce any supplementary oscillation 
phenomena, e.g. relaxation oscillation or self-pulsations originating from a cross modulation 
of the carrier density. The self-pulsation frequency shows up as an enhancement of the 
oscillation of the laser emission if the two frequencies coincide. This allows an accurate 
determination of the SPF frequency close to threshold. 

Figure 2 a shows the optical power and SPF as a function of the bias current. The L-I 
curve has been recorded using an integrating sphere. It is assumed that all emitted power is 
collected. The laser is characterized by a threshold current of ~ 44 mA and a slope efficiency 



of 0.22 mW/mA. The SPF varies from 1 to 4 GHz in a bias current range of 46 to 64 mA, 
which was the maximum injection current we could reach with these devices. In the region of 
the lasing threshold the experimental values present a square-root like behavior dependence 
reminiscent of standard relaxation oscillations as exhibited by a CW-semiconductor laser. 
For bias currents above 55 mA this dependence was no longer observed and the SP behavior 
appears to have a more linear dependence on the bias current. 

III. THEORY 

In this section we use a simple model to explain the observed bias-current dependence 
of the SPF. The investigated laser has a narrow-stripe geometry, which can be modeled in a 
straightforward way using rate-equations || for the optical intensity S (suitably normalized 
to represent the number of photons in the cavity), the number of electron-hole pairs i\q in 
the pumped region, and the number of electron-hole pairs N 2 in the unpumped (absorbing) 
region: 

— = [^(Aq - N tl ) + g 2 (N 2 - N t2 ) - K }S + R sp + F s (t), (la) 

dNl J Nl (M m W *1 ~ vN * nw 

— - = giiNi-NajS , lb 

at e r s T 12 

dN 2 N 2 Ni/v-N, 

—7- = 92{N 2 - N t2 )S + . (lc) 

at t s i 21 

where gi (g 2 ) is the gain coefficient at the transparency number N n (N t2 ) in the pumped 
(unpumped) region, k is the total loss rate. R sp = j3 sp rj sp Ni/T s is the spontaneous emission 
rate, r] sp is the spontaneous quantum efficiency, (3 sp is the spontaneous emission factor and t s 



is the carrier lifetime. Fg(t) is a delta-correlated Langevin noise source JT5| with correlation 
< F s (ti)Fs(t 2 ) >= 2R sp S5(ti — t 2 )/r s , J is the bias pump-current, e is the elementary 
charge, v = V\/V 2 is the volume ratio of pumped and unpumped region, T 12 is the diffusion 
time from the pumped region to the unpumped region, and T 2 \ is the diffusion time from 
unpumped to pumped region. These two diffusion times are interrelated through the volume 
ratio v 1J|: 



,= * = *£. (2) 

v 2 r 21 ^ ' 



Our model (p~a] - p^) is a simplification of the model used in Ref. 0, where the carrier de- 
pendence of the carrier lifetime t s is taken into account using the well-known second order 
expression for r s _1 in the carrier number Nj. Here, we neglect this dependence for the 
moment, as it simplifies the analytical work and qualitatively gives similar results. 

Using the parameter values listed in Table 1, Eqs. ( |Ta] - |Tc| ) are numerically solved with a 
standard algorithm ||. In Figure 3 we show the resulting SPF-J curves, with and without 
spontaneous emission noise. Each value of the curves is calculated from an average over 10 3 
pulses. It is seen that the observed kink in the SPF-J curve is the result of spontaneous 
emission noise. There is a shift of the kink towards larger currents upon increasing the 
spontaneous emission level. For the values of table 1 and (3 sp = 1.3 x 10~ 6 , J xover ~ 82 mA. 
It should be noticed that we do not expect a quantitave agreement between experimental 
and numerical results, since the model neglects important effects, such as gain saturation. 
Nevertheless, the qualitative trends are well reproduced allowing us to physically understand 
the origin of the experimental features. 

Figure 4 shows time traces of the intensity for different bias currents. Clearly, the in- 
terpulse intensity drastically increases with current in the vicinity of the kink. For currents 
J « J X over the interpulse intensity is well dominated by the spontaneous emission (panel 
a)), while for currents J > J xover spontaneous emission does not affect the intensity signifi- 
cantly. The kink-current J XO ver can be defined as the highest current at which the interpulse 
intensity is dominated by spontaneous emission noise. As can be seen in Eq. flla]) sponta- 
neous emission increases the intensity generation rate with an amount R sp . The effect of this 
on the self-pulsation process depends on the generation rate through stimulated emission 
Rstim = [gi(N 1 -N tl )+g 2 (N2-N t 2)]S. For currents J < J xover , R sp > R s tim in the interpulse 
region while for J > J XOV er, the contrary happens. Therefore, the kink pump current J XO ver 
could be mathematically identified through 

J^sp = ■ti' S ti m {J xover ) j \d ) 



where the current dependence of R s um reflects the need to solve Eq. ([|) implicitly using all 
three equations ( |Ia] - |Tc| ) at the time at which the intensity reaches the minimum. 

The long-dashed curve in Fig. 3 is obtained by putting R sp = in Eq. flTa]). In that 
situation, the interpulse intensity becomes extremely small upon decreasing the pump cur- 
rent. The smaller the interpulse intensity becomes, the longer it takes for the absorber to 
reach transparency. When including noise (R sp ^ 0), the interpulse intensity remains at a 
much higher level in the same pump current interval because of the spontaneous emission 
rate R sp . This will significantly increase the speed with which a new pulse is generated after 
the previous one has depleted the absorber. We note that the Langevin noise source Fs(t) 
in Eq. ( |laj ) is responsible for the timing jitter of the pulses. In the region J < J xover a single 
noise event in between pulses may significantly delay or advance the birth of the next pulse, 
causing substantial jitter. For pump currents above the cross-over, the relative effect of 
the noisy events, and hence the jitter, is much smaller. The existence of two pump current 
regions with very different jitter characteristics was also found in Ref. ||. 

In figure 5 the maximum pulse intensity (S max ) and the minimum interpulse intensity 
{Smin) vs. the bias current are shown. An abrupt change (note that the scale in panel c) is 
logarithmic) of S m i n can be seen at J XOV er (while S max takes it maximum value). The kink- 
current J XO ver is therefore identified as the highest current at which the interpulse intensity is 
dominated by spontaneous emission noise. The kink also denotes the boundary between two 
regimes that can be described as follows: For currents larger than J xover the self-pulsation 
has the character of undamped RO, while for currents below this value clear self-Q-switching 
takes place. Obviously, for currents J > J xover the absorber is not depleted deeply enough 
to cause a Q-switch: as soon as transparency is reached, the absorber is bleached but the 
pump is strong enough to prevent total bleaching. For currents J < J xover , the pump is 
small enough to allow total bleaching of the absorbing regions, after which the number of 
electron-hole pairs in the absorbing region has to start all over again. No bifurcation in the 
usual sense can, however, be attributed to this critical current. 

In the next section, we will look at the relationship between ROF and SPF in more 



detail. 



IV. RELAXATION OSCILLATIONS AND SELF-PULSATIONS 

In the previous section we introduced a simple model which provides an explanation for 
the peculiar cross-over region in terms of the average level of spontaneous emission. Here we 
will put our numerical findings in an analytical framework, which leads to a clearer picture 
of the self-pulsation characteristics. 



This is achieved by solving for the CW solutions of Eqs. (|Ta| - |Tq) and investigating their 
stability properties. First we look for laser threshold, which is defined as the circumstance 
for which the trivial solution (S = 0) looses stability in the absence of spontaneous emission. 
We therefore put R sp = in Eq. ([la]) and obtain: 

AT 9i N t i + g 2 N t2 + k , 

Nth ~ g. + V * ( 4 ) 

yi ~ T 12 + VT S 
Jth AT \ l l V Ts 1 (^ 

e t s T 12 T 12 T 12 + vt s 

Using the parameters listed in Table 1, we find J t h = 44.53 mA. 

In total, Eqs. ( |J~a] - |lcD have three possible CW solutions. Below threshold, only the 
solution with S = is physically meaningful (the other two have negative power). At 
threshold the solution S = becomes unstable while one of the other two becomes stable 
with positive power. This is found after performing a standard linear stability analysis, which 
yields for every CW solution a set of (complex) characteristic exponents A = A r -H'Aj. When 
any of these exponents has a positive real part (A r > 0), the CW solution is unstable. The 
imaginary part Aj denotes the frequency with which perturbations initially will grow. Figure 
6 shows how the real parts of the characteristic exponents of the relevant CW solution vary 
with bias current. The CW solution is found to be unstable on the interval 44.556 < J < 92 
mA. For bias currents J > 92 mA, stable CW emission is found. On the other side of the 
interval, a more complex behavior is found. At J = 44.53 mA, the CW solution is stable, 
but looses its stability already at J = 44.556 mA. This sequence of bifurcations from the 
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nonlasing (S = 0) state to self pulsation occurs in a very narrow range of currents around 
threshold. Thus the sequence will be experimentally very hard to resolve due to different 
noise sources; the laser will seemingly begin to oscillate as soon as it crosses threshold. 

Thus, our model ( |Ta] - |ic|) shows that there exists a CW solution that looses stability at 
J = 44.556 niA and regains stability at J = 92 mA. In between these values, the CW state is 
unstable, as indicated by a complex conjugate pair of characteristic exponents with positive 
real parts (Hopf- instability). The region of instability coincides obviously with the region of 
self-pulsating behavior, and is bounded by two Hopf-bifurcations. When the laser operates 
at a bias current 44.556 < J < 92 mA, small perturbations to the CW state in question 
initially grow as exp [(A r + i\i)t], i.e., with angular frequency Aj. The linear stability analysis 
does not provide any information on how this initial growth will saturate. Numerical results 
from Eqs. ( |I~al - |lcl) show that the resulting SPF is always smaller than the ROF Aj/27r. This 
is illustrated in Fig. 3. Both frequencies meet at J = 44.556 mA and at J = 92 mA, the 
two Hopf bifurcation points. In the former case it means that the SPF must increase when 
coming from higher bias currents to reach the RO value. However, this increase only occurs 
in a very small range of currents so that it would be very hard to observe in the experiment. 

It should be noted that a different scenario is found in Ref. |5|]. There, the carrier 
lifetime t s is considered to be carrier dependent, to account for the radiative, non-radiative, 



and Auger processes p2| : 



T-j(N j )=A nrJ + B j N j + C J Nf, (6) 

where j '• = 1 denotes the pumped region and j '• = 2 denotes the unpumped region. This carrier 
dependence is considered necessary because during the strong pulsations, large variations in 



the carrier numbers N~ may occur [22 



It was found in Ref. JJ that the carrier dependence of t s> j(Nj) plays a significant role 
around threshold. This is in sharp contrast with the well-known CW edge-emitting lasers 
where N is clamped immediately above threshold. The kink region, lying far above thresh- 
old, is not affected significantly by taking into account the carrier dependence of r s . This 



illustrates the robustness of the cross-over behavior. At the high end of the self-pulsation 
interval, also a Hopf bifurcation is found, but the dynamics at the low end differs from the 
one discussed here. First of all, there is no window of stability just after threshold. Figure 
7 shows the location of the various CW solutions as a function of bias pump current. The 
S = solution (horizontal solid line) is only shown for currents where it is stable. It looses 
stability at J = J t h- Around J = 0.85Jth two CW solutions are born out of a bifurcation. 
Both CW solutions are linearly unstable. This is in contrast with the model discussed above 
where the upper branch CW solution is stable in a short pump interval after its birth. The 
bifurcation which starts the self-pulsation is not a Hopf one but a homoclinic bifurcation 
(collision of a limit cycle and saddle). Thus, in the model of Ref. || self-pulsation occurs in 
a region bounded by a Hopf-bifurcation on the high bias side and a homoclinic bifurcation at 
the low bias side. This type is not uncommon in (passive Q-switching) self-pulsating lasers 
with saturable absorbers JT^-fJl|]. 



V. CONCLUSIONS 

We have investigated, both experimentally and theoretically, the dependence of the self- 
pulsation frequency of semiconductor CD lasers upon changes in the bias current. A distinct 
kink is found in this dependence, which is investigated using a rate-equation model. We 
have identified that the kink is caused by spontaneous emission, whose average intensity sets 
a lower bound on the emitted laser intensity and thereby on the average intensity, which 
determines the relaxation oscillation frequency. 

From our analysis we conclude that below the crossover point, the self-pulsations be- 
have as passive Q-switching oscillations while above the crossover the behavior approaches 
undamped relaxation oscillations. 

The relationship between the relaxation oscillation frequency and the self-pulsation fre- 
quency is investigated by means of a small signal analysis. We observe that the relaxation 
oscillation frequency so obtained is an upper limit for the self-pulsating frequency. It is also 
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found that self-pulsation occurs in a bias current interval bounded by two Hopf bifurcations. 
A small window of stable CW emission is found very close to the laser threshold in the ab- 
sence of spontaneous emission. The model of Ref. [||] does not show such a window of stable 
emission terminated by a Hopf bifurcation, but a homoclinic bifurcation is responsible for 
the onset of the self-pulsating behavior. However, for the lasers we used in the experiment, 
such differences between the models are irrelevant since they occur in a very small range of 
currents too close to threshold to be resolved. These results raise an interesting question on 
the nature of the bifurcation at the lower side of the self-pulsation interval. 

Note added: A bifurcation analysis of the Yamada model neglecting interstripe diffusion 
has been published |[23| . The bifurcation scenario is different from the one for our model 
(P~a])-(pr^) (and closer to the one in ||). However, these differences in the deterministic 
behavior have no physical relevance, since they appear in parameter domains for which the 
dynamics is dominated by noise. 
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TABLES 



TABLE I. Meanings and values of the parameters in the rate equations 



Parameter Meaning 



Value 



Units 



gi Gain parameter of active region 

<72 Gain parameter of absorbing region 

k Inverse photon lifetime 

t s Carrier lifetime 

f3 sp Spontaneous emission coefficient 

rj sp Spontaneous quantum efficiency 

Nti Carrier number at transparency (active region) 

Nf2 Carrier number at transparency (absorbing region) 

Jth Threshold current of the solitary laser 

a Linewidth enhancement factor 

v Ratio of the active and absorbing volumes 

T\2 Diffusion time 



4.7 x 1(T 9 


ps -1 


1.5 x irr 8 


ps -1 


0.4 


ps -1 


1.1 


ns 


variable 


dimensionless 


0.33 


dimensionless 


6. 10 7 


dimensionless 


1.2 10 8 


dimensionless 


44.53 


mA 


5 


dimensionless 


0.115 


dimensionless 


2.1 


ns 



15 



FIGURES 

FIG. 1. Experimental set-up. 

FIG. 2. Experimentally observed bias current dependence of the emitted laser power (solid 
line) and the self-pulsating frequency (plus signs). 

FIG. 3. Self-Pulsating Frequency (SPF) vs. Bias Current, obtained by numerically solving 
Eqs. (0-IT3). The solid line indicates the value of the relaxation oscillation obtained from the 
small signal analysis. Long-dashed line: SPF in the absence of noise; dash-dotted line: SPF with 
[3 sp = 1.3 x 10~ 4 ; short-dashed line: SPF with fi sp = 1.3 x 10~ 6 ; and dashed-triple-dotted line: 
P sp = 1.3 x 1(T 8 . 

FIG. 4. Time traces of the intensity as a function of the bias current: a) J = 45 mA, b) J = 80 
mA, c) J = 84 mA. (3 sp = 1.3 x 1(T 6 

FIG. 5. Self pulsating frequency (top a)), maximum pulse intensity S max (middle b))and min- 
imum interpulse intensity S m i n (bottom c)) vs. bias current for (3 sp = 1.3 x 1CP 6 . 

FIG. 6. The real part A r of the larger characteristic exponent of the CW solution with positive 
intensity as a function of bias current. The inset, where the real part of all eigenvalues is included, 
shows the tiny window of stable CW emission just after threshold. 

FIG. 7. The intensity of the three CW solutions as a function of bias current, when the carrier 
dependence of the carrier lifetime is included. 
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